Predicted electric field near small superconducting ellipsoids 
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We predict the existence of large electric fields near the surface of superconducting bodies of ellip- 
soidal shape of dimensions comparable to the penetration depth. The electric field is quadrupolar in 
nature with significant corrections from higher order multipoles. Prolate (oblate) superconducting 
ellipsoids are predicted to exhibit fields consistent with negative (positive) quadrupole moments, 
reflecting the fundamental charge asymmetry of matter. 

PACS numbers; 



Electrostatic fields cannot exist inside normal metals 
because the mobile electrons will move along any exist- 
ing electric field lines so as to minimize their potential 
energy, resulting in an equilibrium configuration with no 
electric field. Any kinetic energy gained by the electron 
in this process will be dissipated through inelastic scat- 
tering mechanisms. 

Electrons in superconductors are also mobile (even 
more so than in normal metals) , hence one is accustumed 
to think that no electrostatic field can exist inside super- 
conductors cither. However, we suggest that such pre- 
conception may actually be incorrect for the following 
reason: in a superconductor there are no inelastic scat- 
tering mechanisms that will dissipate the kinetic energy 
of a superfluid electron moving along an electric field line. 
As a consequence, as the electron moves along the field 
line and decreases its potential energy it will gain kinetic 
energy; when it reaches the minimum potential energy 
its kinetic energy will be maximum so it will 'overshoot' 
and move back to a region of higher potential energy, 
leaving the electric field unscreened. Thus we argue that 
at least in principle the existence of electrostatic fields 
inside superconductors is not ruled out by basic physics 
principles. 

On the contrary, we have recently proposed that an 
electrostatic field indeed exists inside superconductors 
0, 13 • The proposal is based on the theory of hole 
superconductivity 4] , that predicts that negative charge 
is expelled from the interior of the superconductor to- 
wards the surface, resulting in a net positive charge den- 
sity in the interior, a negative charge density near the 
surface, and an outward pointing electric field inside the 
superconductor. 

In fact, the possibility of an electrostatic field in a su- 
perconductor follows directly from London's equation for 
the supercurrent J in the presence of a magnetic vector 
potential A|5j 



J = -^1 

TOeC 



(1) 



and Faraday's law in the form that relates the electric 
field E to the electric potential (j) and the magnetic vector 



potential: 



E = -V(/)- 



IdA 
c~dt' 



(2) 



Taking the time derivative of Eq. (1) and using Eq. (2) 
yields 



at rtip 



(3) 



indicating that an electric field that derives from a poten- 
tial will not lead to a time variation of the supercurrent, 
and in particular will not generate a supercurrent if one 
is not present initially. 

The electric potential in the interior of the supercon- 
ductor is assumed to obey the fundamental equation^ 



(f) = -AnXlp{f) + Mr) 



(4) 



which is derived from the London equation under the as- 
sumption that the magnetic vector potential obeys the 
Lorenz gauge condiction^, Q . We postulate Eq. (4) to 
describe the electric potential (j){f) and charge distribu- 
tion p(f) in the interior of all superconductors, with Al 
the London penetration depth and 4>oi'i^ the potential 
originating from a uniform positive charge density po'. 
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where the integral is over the volume of the supercon- 
ducting body. The charge density po is a function of 
parameters describing the superconducting material, of 
the dimensions of the body, and of temperature It 
originates in the 'undressing' of carriers as the system 
goes superconducting, as described in the theory of hole 
superconductivity In particular, it increases as the 
temperature is lowered below Tc and the superfluid den- 
sity increases, and it is larger for superconductors with 
higher Tc- 

The electric potential also obeys the usual Poisson 
equation 



-47rp(r) 



(6) 
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in the interior of the superconductor, and Laplace's equa- 
tion in the exterior. Furthermore we assume that (jj^f) 
as well as its normal derivative are continuous across the 
surface of the superconductor, i.e. that no surface charge 
density exists. 

Equation (4) predicts that p(r) = po deep in the inte- 
rior of superconductors of dimensions much larger than 
the penetration depth, and an excess of negative charge 
within a layer of thickness \l of the surfacel^ . The solu- 
tion of Eqs. (4)-(6) was discussed in Ref.Q for the case 
of spherical superconductors. For that case, no electric 
field exists in the exterior of the superconductor and the 
proposed scenario cannot be detected. Fortunately, the 
situation is different for superconductors of non-spherical 
shape. 

We consider here ellipsoids of revolution, defined by 



z 
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= 1 



(7) 



with p, z cylindrical coordinates. The dipole moment of 
any charge distribution with the symmetry of the ellip- 
soid is zero, but the quadrupole moment 



(8) 



is not. For the uniformly charged ellipsoid with charge 
density po the quadrupole moment is 




FIG. 1: Electric field outside a prolate ellipsoid of dimensions 
a = \, h = 1.5. London penetration depth is Ai = 0.5, and 
Po = 1- The arrows point in the direction of the electric field 
and their length is proportional to the magnitude of the field. 
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so that it is positive (negative) for prolate (6 > a) (oblate 
{b < a)) ellipsoids for positive po. 

We solve the differential equation resulting from Eqs. 
(4)-(6) 



(f) = Aiv2</)(rO + 0o(^ 



(10) 



numerically in the interior of the ellipsoid using the 
GENCOL algorithmlij. Initially we assume an arbitrary 
boundary condition consistent with overall charge neu- 
trality, for example that the normal derivative d(j)/dn = 
everywhere on the surface (Neumann boundary condi- 
tion), which implies non-existence of electric fields in the 
exterior of the superconductor. However d4>/dn = on 
the surface implies ((> is constant on the surface. The first 
iteration of the procedure yields a on the surface that 
is not constant (except for the case b — a), implying that 
the original assumption of no electric field in the exterior 
was incorrect. 

To achieve self-consistency we obtain the electric po- 
tential in the exterior of the superconductor using 



J3J 



p{f 



(11) 



where p{r) is obtained from Eq. (4) using the solution 
obtained for (j)(f) in the interior. For the next iteration 
step we can use as boundary condition for the interior 



problem the normal derivative d4>/dn obtained from the 
exterior problem or the potential itself, obtained from Eq. 
(11) on the surface. After a few iterations the procedure 
converges to a unique solution (/)(r) which is continuous 
and has continuous derivatives across the surface. 

For a prolate ellipsoid with b/a — 1.5 and London pen- 
etration depth Xl/ci =^ 0.5 the quadrupole moment of 
the resulting charge distribution is Q = —.485, in units 
so that pq = a = 1. For comparison the quadrupole mo- 
ment of the bare positive ellipsoid for this case is Qo = tt . 
Figure 1 shows the electric field configuration outside the 
ellipsoid obtained (the electric field inside is not shown 
because it is much larger). The electric field points out 
near the surface for p ~ a, z ~ and points in near the 
surface for p ~ 0, z ~ 6. For an oblate ellipsoid the situa- 
tion would be reversed. The tangential component of the 
field on the surface shown in Fig. 1 exerts a force on the 
electrons in the direction of smaller \z\. As argued earlier, 
we interpret that electrons do not rearrange according to 
this force to screen the field because they will increase 
their kinetic energy when they move in the direction of 
the force and decrease it when they move opposite to it. 
The magnitude of the electric field decays rapidly away 
from the surface since it is of quadrupolar nature. 

Figure 2 shows the obtained charge density p(r) in the 
interior along the axis z = and p = and the electric 
field in the interior along these axis for the parameters of 
Fig. 1. Because Al is comparable to the size of the body. 
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FIG. 2: For the parameters of Fig. 1: (a) Charge density in 
the interior of the superconductor along the horizontal axis 
plotted versus p/a (curve labeled pja) and along the vertical 
axis plotted versus zjh (curve labeled z/b) and on the bound- 
ary plotted versus 6/{tv/2), with 9 = tan~^ {z/b)/{p/a)); note 
that the negative charge density near the surface is larger in 
magnitude along the z direction, (b) Electric fields in the in- 
terior along the p/a (dashed) and z/b (dot-dashed) directions. 
Note that the electric field along the z direction changes sign 
near the surface, and that the electric fields are finite at the 
surface. 



p{r) does not reach the value po {—1) deep in the interior. 
Contrary to the case of the sphere discussed in ref. 0, 
here the charge density and electric fields are different in 
the different directions. Furthermore the negative charge 
near the surface cannot fully screen the electric field as in 
the case of the sphere, so that the field is finite at the sur- 
face of the ellipsoid and leaks out. Note that the obtained 
negative charge density near the surface is larger in the 
direction of the larger axis of the ellipsoid. This is found 
to be the case generally for both prolate and oblate ellip- 
soids, which implies that the quadrupole moment will al- 
ways be negative (positive) for prolate (oblate) ellipsoids. 
This would even be the case if the charge density near 
the surface was equal in the directions of large and small 
axis, but this feature enhances the effect even more. For 
smaller values of the penetration depth the charge den- 
sity is larger in the interior and more negative near the 
surface as in the case of the sphereQ, and a finite electric 
field always occurs at the surface except for a = b. 

Figure 3 shows the electric potential as well as the nor- 
mal and tangential components of the electric field on the 
surface of a prolate ellipsoid for the parameters of Fig. 1. 
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FIG. 3: Electric potential (a) and electric field in the normal 
(b) and in the tangential (c) directions on the surface of the 
ellipsoid for the parameters of Fig. 1 plotted versus 6/{n/2), 
with 9 = tan~^ {z/b)/{p/a)). The solid and dashed-dotted 
lines give the results obtained from the inside and the outside 
solutions, as described in the text; the dashed lines give the 
corresponding potential and fields for a quadrupole at the 
origin of the same magnitude as the quadrupole moment of 
the charge distribution for this case, Q = —0.485. 



The agreement of the results from the inside and outside 
solutions indicates that the numerical iteration procedure 
has converged. The electric potential is larger near the 
equator {6 — 0) than near the poles {6 — tt/2), consis- 
tent with a tangential electric field in the polar direction 
and a normal field that points out near the equator and 
in near the poles. The results differ substantially from 
what would be obtained from a pure quadrupole with ro- 
tational symmetry along the z axis, where the potential 
is given by 



(12) 



These results are shown in Fig. 3 as dashed lines for com- 
parison. This implies that there is substantial contribu- 
tion from higher order multipoles. Moving away from the 
surface we find that the pure quadrupolar behavior pre- 
dicted by Eq. (12) rapidly becomes dominant. It is easy 
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to obtain the contribution from the different higher order 
multipoles by either comparing the potential outside for 
different distances or by computing integrals of the ob- 
tained charge distribution inside. A precise experimental 
study of the field distribution in the neighborhood of su- 
perconducting ellipsoids should give detailed information 
on the charge distribution inside the superconductor. 

As Al decreases compared with the dimensions of the 
body the electric fields increase, and the electric field in 
the interior approaches the one of a uniformly charged 
ellipsoid of charge density po, given by Ep = Cpp^p, 
Ez — CzPqz, with Cp,Cz scale-independent geometrical 
constants. The electric fields at the surface are a fraction 
(whose value depends on b/a) of the fields at the sur- 
face of the uniformly charged ellipsoid. For example for 
b/a =1.5 the limiting values for the fields at the surface 

are approximately E'^ - 0.13£;;"''^, E^ 0.3^;™"^, 

with E'^'^^ , E"^""^ the maximum values of the fields in- 
side, i.e. CpPoa and CzPob respectively. The difference in 
potential at the surface between the poles and the equa- 
tor reaches a limiting value which for the case b/a = 1.5 
is approximately A(j) ~ 1.2Epa. 

As discussed in ref. 0, energy considerations show 
that Po decreases inversely with the linear dimension of 
the sample. Consequently the electric fields at the surface 
first increase as the sample size increases and then reach 
their limiting values once the sample is much larger than 
the penetration depth. The potential difference A(/) at 
the surface increases linearly with the linear size of the 
sample when the fields have reached their limiting values. 
We estimated in ref. Q the maximum electric field in the 
interior of both high Tc cuprates and Nb to be of order 
lO^V/ cm, which according to the above discussion then 
yields for ellipsoids with b/a ^ 1.5 electric fields near 
the surface of order lO^V/cm for samples larger than the 
penetration depth, and potential difference between polar 
and equatorial points on the surface A0 ^ lif'V x a{cm), 
with a the linear dimension of the body. 

Such large electric fields are not seen near the surface 
of macroscopic superconductors. We propose that this is 
so because when the potential difference A0 on the sur- 
face becomes larger than the work function W , electrons 
will 'pop out' of the superconductor near the region of 
low potential and migrate to the region of high potential 
on the surface, and the resulting electronic layer outside 
the superconductor will screen the electric field so that it 
becomes unobservable. For the example discussed above 
with a work function W ~ ^eV this implies that the 
predicted large electric fields will only be observed for 
samples with linear dimensions smaller than a ^ 5000A. 
We emphasize however that this estimate is very rough 
and depends on intrinsic parameters of the superconduct- 
ing materialQ and on the shape parameters. The maxi- 
mum sample size for which the fields are unscreened will 



increase as the temperature increases (since po will de- 
crease) and the magnitude of the fields will decrease, such 
that the potential difference on the surface never exceeds 
the work function. 




FIG. 4: Superconducting ellipsoids will tend to bunch into 
larger spherical shapes. The arrows indicate the direction of 
the electric fields at the surface. 



Note also that an attractive or repulsive electric force 
will exist between small superconducting ellipsoids de- 
pending on their relative orientation. A set of small 
superconducting ellipsoids in close proximity will lower 
their energy by adopting the proper orientation and com- 
ing close to each other, thus bunching into spherical ag- 
gregates as shown in Fig. 4. We suggest that this may 
be related to the remarkable experimental finding of for- 
mation of superconducting balls reported by Tao and 
coworkers 0. 

Unlike the conventional description that is non-Lorentz 
invariant, Eq. (1) with A in the Lorenz gauge leads to 
a relativistically covariant description of the electrody- 
namics of superconductors 0, IT(i lllj |. which we do not 
believe violates any basic physical law nor contradicts 
any known experimental fact. The resulting electric field 
in the interior of superconductors does not lead to a time 
dependent supercurrent because Eq. (3) replaces the con- 
ventional relation between current and electric field for 
'perfect conductors' usually assumed to be valid for su- 
perconductors. Further experimental consequences and 
the relation with microscopic theory will be the subject 
of future work. 
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